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Abstract 

We present a sufficient condition for the presence of spontaneous magnetization for the Ising 
model on a general graph, related to its long-range topology. Applying this condition we are able 
to prove the existence of a phase transition at temperature T > on a wide class of general 
networks. The possibility of further extensions of our results is discussed. 
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INTRODUCTION 

Since the original paper by Isingjlj], in which it was proved that the Ising model on an 
infinite linear chain would not show a phase transition, a huge amount of research has been 
conducted on the subject. The first phase of this fruitful line investigated regular systems, 
and after the articles by PeierlsjsJ and Onsager|3] it became clear that regular lattices in d 
dimensions would magnetize when d > 2. 

In a second phase, a large number of fractals was investigated^ -8], mainly via the so- 
called renormalization group techniques, to discover that, although no rigorous theorem 
has been proven, those (and only those) fractals which have an infinite minimum order of 
ramification display spontaneous magnetization. 

In the same years, fundamental analytical results were obtained for disordered structures 
embedded in Euclidean lattices, applying percolation theory concepts 



More general graphs 



10Hl4j have become increasingly popular in the last twenty years: 



the main difference with the previous cases is that the metric structure of the embedding 
space ceases to play an essential role, as in general a graph is a topological structure which 
is not necessarily embeddable in a finite dimensional Euclidean space. The absence of 
translational invariance and scale invariance makes general graphs very difficult to study, as 
ad hoc techniques must be employed, that usually admit no straightforward generalization. 

An important result would be the identification of a simple parameter, capable of deter- 
mining whether the Ising model on a given graph exhibits a phase transition: we present 
here a theorem stating a sufficient condition for a graph to exhibit spontaneous magnetiza- 

n n 

tion, which is a generalization of the classic Peierls-Griffiths theorem[2|, [15[ for the square 
lattice. While interesting works, employing the same basic techniques as the Griffiths theo- 
rem, have been proposed for higher-dimensional lattices, typically stemming from the paper 



of Dobrushin 
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17], such as the profound contribution by Isakov|l8l| and the extensions to 



non-symmetric situations treated in Pirogov-Sinai theory|19|. or directly from the paper by 
Griffiths, as in Lebowitz and Mazel 20J , nothing applying to inhomogeneous networks and 
arbitrary graphs has yet emerged, and our contribute aims essentially at filling this gap. 

The reason why the modulus of magnetization (|M|) is considered is that it is indissolubly 
tied to the long range order of the graph: it is easy to prove that, when the external field is 
zero, stating (\M\) > e > is equivalent to the existence of a non-zero measure subset of 
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all the correlation functions such that all of its members are greater than a small constant 
e' > 0. 

In the following, we first present the concepts of open and closed borders in a graph for 
later use; we then define the ferromagnetic Ising model on a general graph and derive the 
equivalence of the sum over configurations and the sum over different borders. Next we prove 
a theorem stating a sufficient condition for a graph to exhibit spontaneous magnetization. 
Because of the technical nature of the theorem, we thoroughly examine its more and less 
immediate consequences for a wide range of different graphs. Lastly, we discuss our results 
and the current comprehension of the mechanism of spontaneous magnetization on graphs 
for the Ising model. 



OPEN AND CLOSED BORDERS IN A GRAPH 



A graph Q is a pair (V, £), where V is a countable collection of vertices and £ C V x V 
is a set of unoriented bonds between points. Any pair Q' = (V',£'), such that V' C V, 
£' contains only links between elements of V and £' C £, is called a subgraph, and it's 
denoted Q' C Q. We will restrict our attention to those graphs whose coordination number 
Zi, representing the number of bonds in £ having one extremum in i, is uniformly limited: 
an integer zuax > exists such that Zi < zuax for all i 6 V. 

We now define a path 7, between two points i and j, as a collection of consecutive bonds 
of £, where consecutive means that each pair shares a vertex with the next one: 

7 = {(ih), (hh), (h-i,j)} ■ 

Directly associated to the concept of path, the chemical distance between two points % 
and j is defined as the length of the shortest path connecting them. The chemical distance 
straightforwardly induces the so-called intrinsic metric of the graph. 

The intrinsic fractal dimension df rac of a graph, defined as the minimum d such that 



A,., the maximum number of vertices included in a Van Hove sphere 21[ of radius r (i.e. 
the set of points within a chemical distance from a given point of no more than r bonds), 
satisfies N r < r d as r — > 00. It differs from the usual fractal dimension in that it refers to 
the topological nature of the graph (i.e. on its natural - chemical - distance), and not on 
the metric structure of the space into which the graph is embedded. 
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To proceed we declare what will be considered a border from now on. 

Definition. Given a connected graph Q = (V,£), we can define a border B as a set of 
bonds that separates exactly two connected subgraphs. It means that two sets V\^Vi C V 
exist such that 

• V x n V 2 = and V 1 UV 2 = V, 

• any path on Q from a point of V\ to a point of P 2 must contain at least one bond of 
B, 

• a path exists between any two points in V% (i = 1, 2) that doesn't contain any bond of 
B. 

It is noteworthy that the union of two disjoint borders is not a border itself under this 
definition, as it divides the graph into three subgraphs. This is a feature we'll later need to 
avoid overcounting different configurations. 

The intuitive idea of open and closed border is actually an artifact created by our visualiz- 
ing regular lattices as immersed in a finite dimensional real space: the seeming adjacency of 
the vertices creates a contour of the graph, which we use to define closed and open borders. 
The fact is that this contour is heavily dependent on what particular immersion we employ, 
and ceases to exist when we consider the graph for itself. The border in itself has no geom- 
etry whatsoever, since it is just a collection of links, and even the notion of "continuous" 
border, without further specifications, makes no sense from a graph-theoretic point of view: 
in a general graph, a border is just a collection of links that splits it into two parts. We now 
define open and closed borders with respect to an external set of points, as it will be useful 
later. 

Definition. Given a border B and a set of points £ C V , we say that B is closed with 
respect to the external points set £ if either V\ fl £ = or V2 H £ = 0, otherwise B is open. 

For any finite subgraph Q N of a given graph Q, we choose the natural set of external 
points £: 

£ = {i e Q N : E C for some j &G\ Q N } ■ 
Now, given a border B{ that divides Q N into two subgraphs Ai and Cj, we define A4 as 
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• internal if 

— Bi is closed and A; L n £ = 0, or 

— £>j is open and contains fewer elements than Cj, or 

— is open, Aj has the same size as Cj and the points in A; L linked to Bi have 
negative spin; 

• external if 

— Bi is closed and Ai fl £ 7^ 0, or 

— £>j is open and Ai has more elements than Ci, or 

— Bi is open, Aj has the same size as Ci and the points in A4 linked to Bi have 
positive spin. 

The reason why we had to select a finite subgraph Qn is that we need to be able to count 
the number of spins in the graph for the previous definitions to make sense. 

THE FERROMAGNETIC ISING MODEL ON A GRAPH 

Let now ^ = ±1 be a spin variable for each vertex % e V. We define the Ising Hamil- 
tonian on a graph as 



where the couplings Jij = Jy L must satisfy < Jy < J Max < 00 for some J Max, and Jij > 
if and only if € C In the following we will set the external field to zero everywhere 



Now that we have presented the terminology we'll be using, we are going to study the 
equilibrium statistical mechanics of the Ising model at inverse temperature j3, and in par- 
ticular the modulus of the magnetization 




(1) 



(hi = 0). 




(2) 



where Z = J2s a .\ e ^tt^d) [ s the partition function, and \V\ is the cardinality of V. 
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Since A\M\ = (|M| 2 ) - (\M\)(\M\) is a variance, 

(M 2 > > ((|M|)) 2 , 

so stating (|M|) = e > implies (M 2 ) > e 2 > 0. On the other hand, since M 2 < \M\, the 
converse is true, so (M 2 ) > and (|M|) > are equivalent. 

As we have now defined the main quantities we'll be studying, our next step is to prove 
that we can substitute the sum over configurations of the graph with a sum over possible 
border classes, that we now define. 

EQUIVALENCE BETWEEN SETS OF BORDERS AND SPIN CONFIGURA- 
TIONS 

Definition. A border class is a class C = {C 1 } of border sets C l = {B{, B l 2 , . . . , B^} , 
where i = 1, . . . , N c , such that 

• B l u n B l v = for all i = 1, . . . , N c and u, v = 1, . . . , N i} 
. u i=i,...,jv i _ B i = u m=i,...,j\r f _ B i i f or all i ; j = 1, . . . , Nc- 

Theorem. To any given border class corresponds one and only one configuration of spins 
on V , once we set the value of a single spin. 

Proof. To prove that, for any border class and a given spin p e V, we can construct a single 
spin configuration, we first choose an arbitrary representative C l = {B\, . . . , B l N } of C and 
set all the spins to the value of p, then for each B\ e C l we flip all the spins of the subgraph 
which doesn't contain p. The result is independent of the order in which we choose the B l k , 
since each spins changes sign once for every border that separates it from the fixed spin p, 
and is independent of the specific i. 

To prove that for any given spin configuration we can create a single border class, we proceed 
as follows: let be the sets of all plus (minus) spins, 

& = {i e G : <7i = ±1} ; 

We now choose the subsets Rf of R ± , so that each Rf is connected, while for all i ^ j Rf 
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and Rf are disconnected; moreover we require that 

= Rf U Rf U • • • U Kjj, 

RfnRf = ® v^, 

We are selecting individual clusters of homogeneous spins, so satisfying the above requisites 
is always possible. Setting now 

dRf = {(a,b)eC:aeRf,b^Rf}, 

the sets B ± = dRf U • • • U dRf are a collection of links each defined unambiguously, and 
furthermore B + = B~ . It may happen that for some i the subgraph Q \ Rf is made of two 
disconnected subgraphs (e.g. when a ring of plus spins is surrounded by minus spins); as a 
consequence dRf is not a border according to our definition. In that case it is possible to 
split dRf into subsets, so that each of them divides Q into two connected subgraphs. After 
dealing in this way whith all the Rf, we are left with a collection of well-defined borders 
dTf, with j — 1, . . . , U + , and dT^ , with k — 1, . . . , LL where [/+,£/_ > S. It is still possible 
that some of the borders Tf~ , while defining exactly the same zones, have no correspective 
in T~ but, since they nevertheless verify B + = B~ , they belong to the same border class, 
completing the proof. □ 

The main consequence of this result is that we can substitute a sum over border classes 
for a sum over configurations whenever needed, and we can infer from the structure of 
the borders some limiting properties for the spins distributions, as we'll see soon. It is 
worthwhile to explicitly notice that, when we pass from a sum over configurations to one 
over borders, and not border classes, we overcount some borders, as there are more than one 
representative of each border class: this is not going to be a problem in the use we'll make 
of this result. 

GENERALIZED PEIERLS-GRIFFITHS' THEOREM 

We can divide the set of all configurations on Q N into two classes: 

• all the negative spins are internal to some border (class Af), 

• at least a negative spin exists that is external to all borders (class V). 
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The second case implies that every positive spin lies inside some border, since it must lie on 
the opposite side of the negative spin which is always external. 

We now restrict our attention to the configurations belonging to the first class, denoting 
by a subscript M the quantities that pertain to it; we can obtain a good estimate of the 
number of negative spins, (N_)jy, as follows: the sign of a spin p is negative if it is contained 
inside an odd number of borders, positive otherwise; we obtain a very naive, yet effective, 
approximation if we consider any spin contained inside at least one border as negative: 
letting P be 1 if p is inside at least a border, otherwise, we can write 

(n-u <E< /P >- 

peg 

We are now to give a reasonable estimate of (P) : take all the configurations C with at least 
one border containing p, call b min the length of the shortest border in C containing p and let 
k be the number of borders containing p, so as to write 

(P) = z^ E e-e* = E E E e ~ m 

C\p inside b min >l fc>l Q\ h rnin 

i k borders 

Now fix bmin and consider the configurations containing k borders: if we remove the shortest 
border from such a configuration C, we obtain a new configuration C with (k — 1) borders 
containing p, each of them at least b m i n long. C will be present in the partition function Z, 
but different configurations with k borders C may give the same C; defining now fi p (b) as the 
number of possible borders of b links containing p, the degeneration induced by removing 
the shortest border is not greater than [i p (b min ). The energy of a configuration C and the 
corresponding Boltzmann factor obey 

if we now limit the sum in the partition functions to those configurations obtained by 
removing a border from the numerator, we can write 

Z- 1 E E ^ ^ /A&mm) e~ 2/3J — b — . 

^ — 1 C\k borders 

In this way the average number of minus spins is bounded by a function depending only on 
the number of borders encircling a given spin: 

(N-)« < E E /A&mm ) e" 2/3J — 6 — ; 
8 



this result states that, no matter what the maximum number of spins you can isolate inside 
a border is, as long as fj> p (b) grows at most exponentially the value of (N-)m can be limited 
at low enough temperatures. An analogous result holds for configurations of class V when 
exchanging the roles of positive and negative spins: 

(N + ) v < E E ^ h ^n) e~ WJmm bmm ■ 

P&G b mi „ 

Let now //(&)) = sup p /i p (6). We can now prove the following theorem: 

Theorem. // on an infinite graph /i(b) < A b , definitely for b > b and for some for some 
A > 0, then the graph exhibits spontaneous magnetization at large enough (3 (low enough 
temperatures). 

Proof. The average modulus of magnetization is \M\ = N' 1 (N + — N_); writing the Boltz- 
mann factor for a configuration C as Pq = -2 -1 e - ^ c \ we can write the following: 

(|M|) = \ M <\ p c = E M C P C - M c p c 

CgvuM ce(vuM)+ Ce(VuM)~ 

= I M cPc - E M c p c) + ( E M c p c - E M c p c ) 

= lj2 M cPc-J2 M c p c) + 2 ( E M c p c- E M c p c) 

\ceM cev J \cev+ CeAf- J 

> (j2 M cPc-J2 M c p c) = i-^(E( iV -)^ + E( iV +)^) 

\CeAf CeV J \CeAT CeV J 



^ 1 - m E E e ~ WJmm bmm 



p£V b min 

-2/8(6-l)J„ 

6>1 



> l-Ae- 2 ^J2^b min )e- 



When the sum on the last line converges, the equation tells us that, for large enough j3 (low 
temperatures), (\M\) is greater than a positive constant, so that spontaneous magnetization 
on an infinite graph is achieved, while in general (\M\) is finite for every N, but can tend 
to zero as N — > oo. □ 

The main problem in employing the previous theorem is determining bounds on fi(b). 
The simplest case in which the hypothesis does not hold is a situation in which for some 
finite b the number of borders surrounding a given point is infinite. As a sound check of the 



validity of the theorem, all the weakly separable 22[ graphs, which do not magnetize, belong 
to this category. 

To further our understanding of the result, we need to present a new parameter. Given 
a subgraph A C Q, we define its external boundary dA as the set of points in Q \ A that 
have a bond to a point in A; denoting the number of vertices in A as \A\ we now present the 

d-l 

isoperimetric dimension dj so as the minimum d such that dA > C- \A\ d . The largest 
set of points encompassable with b links is thus smaller than b^°~ x \ since these points are 
connected, b^o- 1 is also the maximum radius of a set including % with a border b, so the set 
of reachable points, V(b) C V, has a cardinality |V(6)| < b ^o- 1 . 



Given a point p and for each border B, consider B, the collection of vertices contributing 
to B which are on the inside of B with respect to p. As the two are in biunivocal relation 
once p is chosen, counting the borders is the same as counting the vertex borders. 

As a consequence of the previous paragraphs, the following holds: 

Proposition. In a graph with isoperimetric dimension di SO > 1, the number of possible 
borders surrounding p is bounded by 

//*(&) < J2 N i( b ) ^ b-^t=^.N Sup (b), 
i 

where the sum is over the points q which can be enclosed in a border of size b, N q (b) is the 
maximum number of vertex borders of length b, containing p, which can be created starting 
from the point q, and N Sup (b) = sup qeV ^ N p (b). 

A border is connected if the corresponding vertex border is a connected set. We will need 
the following proposition regarding N p (b) to obtain a general result. 

Proposition. The number of connected vertex borders starting from a given point p grows 
at most exponentially with b: N p (b) < const- C b . 

Proof. A tree is a graph that has no loops. From each connected subgraph A we can draw a 
number of different spanning trees, i.e. trees having the same set of points V as the original 
A. For any spanning tree we can construct a path visiting all its vertices in no less than 
b — 1 steps, and in no more than 2zM ax b steps. While the former statement is obvious, we 
now prove the latter using the following algorithm: starting from i, choose link and cross it; 
at each vertex on the path, choose a link not yet crossed; if there is no free link, step back 
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FIG. 1. The Sierpinski carpet allows for disconnected borders, but they can be connected with no 
more with I- b links, so the number of possible borders grows no faster than an exponential with b, 
and the graph magnetizes. 



through the link from which the path first arrived at the vertex. With this algorithm, each 
link is crossed no more than twice, so the path is of no more than 2zMaxb', furthermore, all 
the links are crossed, so, since the boundary is connected, all the points are visited. As a 
consequence, all the spanning trees of b points starting from a vertex i can be constructed 
as paths of b — 1, b, . . . , Izuaxb steps. Since each step can be chosen among at most zu ax 
links, the total number of possible spanning trees, starting from i and made of b vertices, is 
less than 

2z Max b+l _ 6-1 

J>-\ _|_ b , _|_ 2z Max b _ ^Max ^Max 

^Max ^Max ' ' ' ' ' ^Max ~ i ) 

ZMax ~ J- 

and the thesis follows: 

N p {b) < const- C b . 

□ 

When a vertex border is made of more disconnected parts, /i(b) grows exponentially if, 
for all borders, it's possible to connect all the parts using no more than /• b vertices, where 
I is a constant of the graph: in fact in this case to each border of length b corresponds one 
connected vertex border of length between b and b- /, so that 

fil-b+l _ fib 

N(b) <C b + C b+1 + ■■■ + C lb = — — — . 

C — 1 
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Noting that di SO > 1 implies that there is no border length b for which /i p (6) is infinite, 
the previous results can be combined to form the following theorem. 

Theorem. For all graphs with isoperimetric dimension di SO > 1 and vertex borders which 
are connectable with no more than I- b vertices, a finite critical /3 C < oo exists such that for 
all (3 > (3 C spontaneous magnetization is achieved. 

To the latter category belong the regular lattices in d > 2 dimensions and crystals with 
any kind of elementary cells; we explicitly note that for an Euclidean lattice in d — 2 
dimensions we recover the result by Griffiths [15]. In addition, each vertex border can be 
connected with no more than I- b vertices in the Sierpinski carpet too, which therefore 



magnetizes, in accord with the existing literature [23] (see FigJT]). 

Consider now the ladders of infinitely growing height (see Figj2]): they are structures 
described, at any offset n on the semi-infinite base line, by a non- decreasing integer function 
h(n); as long as the isoperimetric dimension of the ladder is strictly greater than one (h(ri) > 
A n a + B , for a > 0, A > and B > 0) the previous arguments apply, so the ladder 
magnetizes; on the other hand, when di SO = 1 a little more work is required: the total 
number V(b) of vertices which can be included in at least one border of length b can not 
grow faster than the number of points on the left of the rightmost border of length 6; the 
latter is at at offset n{b) = max{n'|& = h~ 1 (n')}, so that V(b) satisfies V(b) ~ Ya=i M*) f° r 
large b; when 

B + A log i < h(i) < B x + A x i a 
for some Aq,A\ > 0, Bq, B\ > and a > 1, the volume satisfies 

pe b / A 

V(b) < J di {B l + A\i a 



b/A ° 1 A, 

£l e 2(a+l)b/A 



A 2 (a + 1) 

and since the borders are all connected n(b) is exponential and the graph magnetizes. When 
instead lim„_ s>00 = 0, for all e > and n large enough h(n) < elogn holds; as a 
consequence, 
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V(b) > / di h(i) > di = e~ 
Jo Jo 

holds for all e > 0; in this case the sum /i(6)e _2/3 ^ _1 ^ Jmin diverges for all temperatures, 
so the hypotheses of our theorem are not fulfilled. These results are in agreement with a 



result by Chayes and Chayes[24( about more general structures called d- wedges, where it 
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■■■■ rnnrnmr 

a) 
b) 

FIG. 2. Two examples of growing ladder graphs: in (a) all the borders are connected (remember 
that, following the definition we use here, a border divides a graph into exactly two subgraphs, each 
connected), and as a consequence the graph will magnetize. In (b) too the borders are connected, 
and even if the growth in width of the ladder is very slow (logarithmic) the same result holds. 

is proved that h(n) > logra is both a sufficient and a necessary condition for spontaneous 
magnetization. 

DISCUSSION 

To give a more intuitive interpretation of the theorem we proved, we can proceed as 
follows: if the number of borders grows less than exponentially, we can argue that all of 
these borders will contain a number of spins increasing slowly with the length of the border; 
as a consequence, the formation of large clusters of spins in a magnetized graph will be 
energetically unfavoured, so that the latter will result a stable state. On the other hand, if 
fj, l (b) grows very fast with b, we expect that some of the borders will be far from the vertex 
i, so that more and more vertices will be enclosed in short (low b) borders; this in turn 
means that large clusters of spins can be nipped spending a small amount of energy, so that 
a magnetized graph may be unstable with regard to thermal fluctuations. 

The condition of our theorem is a strong one, in that it investigates a global property 
of the graph. For this reason it can not be a necessary condition for achieving spontaneous 
magnetization: if a graph has a part, which has zero measure in the thermodynamic limit, 
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for which the number of boundaries ^P{b) is greater than any exponential (e.g. a semi-infinite 
line connected to a point on a plane), the hypothesis of the theorem is false but the graph 
clS cL whole can still magnetize. 

An important, yet straightforward, observation is that whenever a subgraph of non zero 
measure exists that is magnetizable, all the graph is magnetizable: in fact all the correlation 
functions, as computed on the subgraph, are smaller than or equal to the corresponding 
ones in the complete graph; when, on the other hand, the graph is formed by a collection 
of zero measure, weakly connected, magnetizable subgraphs (e.g. an infinite collection of 
parallel planes, each connected via a single link to the next one), there is no guarantee that 
(|M|)>0. 

Our result about the Ising model on graphs is a further step towards a full comprehension 
of the mechanism of phase transitions on general networks: together with a sufficient condi- 
tion for the lack of spontaneous magnetization 22] , it allows to ascertain the magnetizability 
of a large number of structures with a minimal amount of computation. 

Further steps extending this work should aim at closing the gap between magnetizable and 
non- magnetizable graphs under the (\M\) definition, in order to identify a condition both 
necessary and sufficient for spontaneous magnetization; another direction of development 
could be to treat non symmetric situations, as in Pirogov-Sinai theory. 



[1] E. Ising, Z. Phys. 31, 253 (1925). 

[2] R. Peierls, P. Camb. Philos. Soc. 32, 477 (1936). 

[3] L. Onsager, Phys. Rev. 65, 117 (1944). 

[4] Y. Achiam, Phys. Rev. B 31, 4732 (1985). 

[5] M. A. Bab, G. Fabricius, and E. V. Albano, The J. Chem. Phys. 128, 044911 (2008). 

[6] J. M. Carmona, U. Marini Bettolo Marconi, J. J. Ruiz-Lorenzo, and A. Tarancon, Phys. Rev. 
B 58, 14387 (1998). 

[7] Y. Gefen, B. B. Mandelbrot, and A. Aharony, Phys. Rev. Lett. 45, 855 (1980). 

[8] P. Monceau and P.-Y. Hsiao, Eur. Phys. J. B 32, 81 (2003). 

[9] J. T. Chayes, L. Chayes, and J. Frolich, Commun. Math. Phys. 100, 399 (1985). 

[10] A. Ceresole, M. Rasetti, and R. Zecchina, Riv. Nuovo Cimento 21, 1 (1998). 



14 



[11] S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, Phys. Rev. E 66, 016104 (2002). 
[12] C. P. Herrero, Phys. Rev. E 65, 066110 (2002). 

[13] R. Jullien, K. Penson, and P. Pfeuty, J. Phys. (France) Lett. 40, 237 (1979). 

[14] T. Regge and R. Zecchina, J. Math. Phys. 37, 2796 (1996). 

[15] R. B. Griffiths, Phys. Rev. 136, A437 (1964). 

[16] R. L. Dobrushin, Dokl. Akad. Nauk 160, 1046 (1964). 

[17] R. L. Dobrushin, Sov. Phys. Dokl. 10, 111 (1964). 

[18] S. N. Isakov, Commun. Math. Phys. 95, 427 (1984). 

[19] I. G. Sinai, Theory of phase transitions : rigorous results (Pergamon Press, New York, 1982). 
[20] J. L. Lebowitz and A. E. Mazel, J. Stat. Phys. 90, 1051 (1998). 

[21] R. Burioni, D. Cassi, and A. Vezzani, in Random Walks and Geometry, edited by V. A. 

Kaimanovich (Berlin, de Gruyter, 2001) pp. 35-71. 
[22] R. Campari and D. Cassi (submitted). 
[23] A. Vezzani, J. Phys. A 36, 1593 (2003). 
[24] J. T. Chayes and L. Chayes, J. Phys. A 19, 3033 (1986). 



15 



